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Abstract 

Discrete analogs of the finite and affine Toda field equations are found corresponding 
to the Lie algebras of series Cn and Cn. 



1 Introduction 

o 
o 

Consider the Toda chain with three discrete independent variables u, v, k (see [T], 0): 
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^fuv—fu — fv+f _ 

1 + he^ v ~^ u 

Here / = f(u, v, k) is an unknown function, and the following notations are accepted. The upper 

index is used to indicate shifts with respect to the third variable k so that f 1 = f(u,v, k + 1) 

and f^ 1 = f(u,v,k — 1). The lower index shows shifts of the first and second variables 

fu = f(u + h,v,k), f- u = /(it - h,v,k), f v = f{u,v + h,k), /_„ = f(u,v - h } k) } h is 

the parameter of the grid such that for small values of h one gets f u — f ~ VhD x f, and 

f v — f ~ yhDyf, and f uv — f u — f v + f ~ hD x D y f, where D x and differential operator with 

respect to x and Evidently the continuum limit h — > of the discrete Toda chain (DTC) (JTJ) 

i/"") ■ gives the usual two-dimensional Toda chain [3] 
O 

Si D x D y f = ef 1 -f -ef-r\ (2) 

in 

Chain (DTC) is closely connected with the well known discrete bilinear Hirota-Miwa equation 

C ' ttuv tufv = t v t u ■ (3) 

> 

Here the unknown t = t(u,v,k) depends also upon three independent discrete variables. Miura 
type transformation t l = eH converts the equation © into the equation ((TJ). 

The chain (JTJ admits evidently periodical closure constraint f(u, v, k) = f(u, v, k + N) which 
reduces it to a finite field equation with dynamical variables f(u, v, 1), f(u, v, 2), ... f(u, v, N). 
The other closure is defined by the degeneration points of the chain. For example, the chain 
(0) truncated by the conditions e~^ u ' v ' ^ = and e^ u ' v ' N+1 ^ = is an integrable finite system 
of discrete hyperbolic type equations. But as a rule finite field reductions of chains are not 
exhausted by degenerate and periodic ones. For instance, the chain (J2J admits a large class of 
reductions connected with the simple and affine Lie algebras [3]. Aforementioned degenerate 
and periodic reductions correspond to the Lie algebras of series A n and A n , respectively. The 
problem of looking for discrete analogs of finite reductions of corresponding to the other 
Lie algebras of the finite grows is still open (see, for example, [5J, and also the surveys [Ij and 
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jHj). In [Hj the examples of truncations of $1} connected with special classes of solutions have 
been discussed. We suggest below discrete analogs of the Toda chain of the series Cn, Cn- 

To shorten the formulae introduce notations z(u, v, k) = he f «- f » + 1 and Af = f uv - 
fu~fv + /) then the chain gets the form z 1 = e A ^z. Without loosing the generality one 
can put h — 1. Really, the parameter h can be made equal to one by the following shift 
f(k)^f(k) + k\nh. 

The main result of the paper is given in the following two statements. 

Proposition 1. The discrete Toda field equation 

Ar 1 - z 



e 



z 1 ' 



e A/J =^, 0<j<N-l, (4) 

Af N _ J_ 

Z N 

admits the Lax pair (see ()54|) - (jo3j) below). In the continuum limit h — > the chain (0J) turns 
into the Toda equation of the series Cn'- 

u(-l) = -u(0), 

U X y = e <k+l)-u(k) _ e u(k)-u(k-l) ^ ^ 
e u(N+l) = Q 

Proposition 2. The discrete Toda field equation 

AT 1 _ z 



e 



e 



z 1 ' 

^— u,v 
■3+1 

zi ' 

N-l 



0<j<N-l, (6) 



Z N 



admits the Lax pair (see (J6l"j) - (}62j) . below). In the continuum limit h — - > it turns into the 
Toda equation of the series Cn- 

u(-l) = -u(0), 

U X y = e "(fe+l)-«(fc) _ e u(k)-u(k-l) ^ ^ 

u(N) = -u(N - 1) 



2 Involutions of the associated linear systems 

The chain admits the Lax pair consisting of two linear discrete equations pQ 

ijj u = e fu ~ f i) - i/j 1 , ip v = if; + e A, " / ~V~ 1 - (8) 

Here the indices of the eigenfunction ip = ij)(u, v, k) are for the shifts according to the rule, 
defined above: ip u = ip(u + 1, v, k), ip v — ip(u, v + 1, k), ip 1 = ip(u, v, k + 1), ip' 1 = ip(u, v , k — 1) 
and so on. 
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Exclude from the system of equations (jEJ) all the shifts of the third variable. As a result one 
gets a linear discrete hyperbolic equation 

if)uv - Vu - e fuv ~ fv (tp v - z^j) = 0. (9) 

Below we will need equations dual to (JSJ) and ©• In order to find the dual equations we use 
the discrete symmetries of the chain Evidently the chain is invariant under involution 
u — ► 1 — v, v — ► 1 — u. However the involution changes the Lax pair (|8|L which takes now the 
form 

2/-u = 2/ + e / -"- / ~y- 1 > y- v = e f -*~ f y-y\ (10) 
The hyperbolic equation turns into the equation 

1 ef v ~f 

y U v--yv — (y u -y) = 0. (11) 

The idea to use two (mutually conjugate) Lax pairs when studying the finite Toda chain 
belongs to the classical work by Darboux |3j. 

By analogy with the continuous case pose the question when two hyperbolic equations (jHJ) 
and (fTTJl are related to each other by a multiplicative transform like y = aip (see, also, jH])? Two 
hyperbolic equations are connected by a multiplicative transform if and only if their Laplace 
invariants are the same jH] . Remind that the Laplace invariants of the equation 

aip + bvp u + apv + dip uv = (12) 

are expressed as 

K, = K 2 = (13) 



Denote through Ki^(u,v,k), K2^(u,v,k) Laplace invariants of the equation © and through 
Ki y (u, v, k), K2 y (u, v, k) - Laplace invariants of the equation (jllj) . Compute all these invariants 
and find 1111 

K ^v = 71' K *u = -> K H = — ' K W = Zi- ( 14 ) 

* * u Z v 

Therefore, if the field variables satisfy the constraint 

z(u+ l,v, k — 1) = z(u,v + 1, k + 1) (15) 

for a fixed value k = k , then the Laplace invariants of these equations will satisfy the following 
conditions 

K ly (u + l,v, k - 1) = K^(u,v, k ), K 2y {u + l,v, k - 1) = K 2 ^(u,v, k ), (16) 
and the conditions 

Kty(u,v-l,ko-l) = K ly (u,v,k ), K 2 ^,(u,v - 1, k - 1) = K 2y (u,v, k ). (17) 

Consequently, there are such functions R = R(u,v), S = S(u,v) that the following relations 
take place 

if)(u, v — 1, ko — 1) = Ry(u, v, ko) and y(u + 1, v, ko — 1) = Si/;(u, v, ko) (18) 
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between solutions of the hyperbolic equations (jHJ) and (fTTf . Multipliers R and S are found from 
the following overdetermined system of linear equations 

R u = Rz_ v e f ^-v- fI v , R v = Rz_ u e f - U ~ f - U ' v , (19) 
= •S't^t — , S v = S — : — , (20) 

Z -u,-v Z 1 

consistency of which is guaranteed by the condition (fTB^) . 



3 Lax pair of the semi-infinite chain 

Equations (fTHj) can be referred to as boundary conditions cutting off the linear equations (JSJ), 
(IK)jl reducing them into the half-line k > ko (or the half-line k < ko). Concentrate on this 
statement. Put first for the simplicity k = 0. Substituting (fTH|) into the equations (jSJ) and (fTT)|) 

for k = one gets 

y _ u = y° + X_ u if°_ u , if° v = if) + H v y° v , (21) 

where X = Se f ~ f - \ H = Re f ~ f ^ . 

The following lemma gives the connection between the functions X and H. 

Lemma 1. Solutions of the system of the equations ()19|) . ()20|) can be chosen to satisfy the 
constraint 

HX = ^—±. (22) 

Proof. It follows from the equations ()19|) and (f2*U|) that functions X and H satisfy the 
similar linear difference equations of the first order 

X V = X^ H v = Hz\e r -*- r \ (23) 

g/u — f u 2 ef v ~f 
X u = X , X v = X , (24) 

z l z 

where := f(u + 2,v,k-l). 

It is not difficult to check that if the condition (}2*2*j) holds at a point (u, v) then it also holds 
at any neighbouring point. 

Let us discuss briefly dynamical variables. Evidently shifts of the eigenfunction ip in the 
positive direction: ip u , ip u 2, ip v , ip v 2, ... as well as shifts of the eigenfunction y in the negative 
direction such as y_ u , y_ M 2, y_ v , y_ v 2, ... are locally expressed trough unshifted (dynamical) 
variables if), if) , ip ±2 , ••• and y, y^ 1 , y ±2 , ... . This is not the case for the shifts on the 
opposite directions which really should be considered as nonlocal variables. Actually these 
variables cannot be expressed through a finite number of dynamical ones. For example, to find 
the variable ifi- u it is necessary to solve the difference equation (with respect to the argument 
k) of the form 

ifj l _ u - e f ~ f - u ip^ u = if). 

Enlarge the set of dynamical variables. In addition to the set of dynamical variables on the 
half-line k > 0, consisting of the functions {ip° , y° , if) 1 , y 1 , ...}, introduce two more variables Y 
and \l/ by setting Y = y®, = if>°_ u . 
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Lemma 2. The shifts Y- u , Y_ v , Y u of the variable Y as well as the shifts ^f v , of 
the variable \I/ are linearly expressed through a finite number of the elements of the enlarged 
dynamical set Y, y°, ip 1 , y 1 , ...}. 

Proof. Some of the equations required follows directly from the definition: Y_ v = y°, 
= ^ ■ Let us shift the first equation (j2*Tj) to the right by one with respect to u and to v and 
rewrite it in the form 

Y U = Y — X v ^° v = Y — X v (^° + H V Y) = Y(l - X V H V ) - X v ^°. 

Due to the Lemma 1 the quantity inside the parentheses is equal to -, hence 

z 

Y u = -Y-X v ^. 

z 

Shift the expression obtained to the left by one respect to u and transform it as 

Y_ u = z_ u Y + X„ UiV z_ u m. (25) 

Shift now the second equation of (|2"TJl by one to the left respect u and v then after slight 
simplification one gets 

= — ~ H_ u y°, 



and 



tf„ = z_ u m + H_ u „z^ u Y. 



(26) 



Lemma 2 is proved. 

The commentary to the lemma. Variables ipi can be shifted upward and to the right on the 
(u, v )-plane while the variables yi - downward and to the left. The new variables are special - 
they can be shifted on three directions. 

Summarize the computations above. Introduce some notations. Denote through P and Q 
infinite dimensional vectors-columns such as 



( Y \ 

o 



y 
y 



y 
v •■■ J 



Q 



( * \ 

</>° 

v ... j 



(27) 



i.e. Po = Y, Qo = ^ and Qi = ip l 1 for i > 1 p = y % 1 . Introduce six infinite dimensional 
matrices, four of which are two diagonal 



A 



( z u 

1 

e / -«~ / 1 







f2 __fl 



v ••• 



B 



( 1 

e^-J -1 

e f -v~ fl -1 

t 












f 2 - v -f 2 



\ 



'.28) 
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c 



( z v 
1 

eft-* 





1 







e f -~ fl 1 
and two others are of the form 

/ V y 



v ••• 



D 



( o 





v •■• 



Ju-f 













-1 

f 2_ f 2 



(29) 









H ^ uv z v 
H v 


V 



(30) 
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The last two matrices have all entries to be zero except the following ones a ,o, ai,o> c o,cb c i,o- 
Make up a system of difference equations 

P_ u = AP + aQ, P_ v = BP, (31) 

Q V = CQ + cP, Q u = DQ. (32) 

Proposition 3. The systems (131)1 and ()32)1 are consistent if and only if their coefficients 
satisfy the semi-infinite lattice of equations 

e Arl = (33) 



for < j < oo. 



(34) 



To prove the proposition it is enough to check that the conditions (P_ n )_„ = (P_„)_ u , (Q u ) v 
(P v ) u hold only together with the constraints (}33|) . (jHIJ). It is clear that the system (JHT1). (132 
is the matrix form of the following system of the scalar equations 

ip{ = e fi ~ fJ vp j - 4> j+1 , yL v = e f -^~ p y j - y j+1 , for j > 0, 



ipi = ip j + e 



fi-f- 



y-u = v 3 + e 



2-1 



2/ 



for j > 1, 



F_ u = + X^ u>v z- U ty, ^ v = z_ u ^ + H_ U)V z„ u Y, 
Y_ v = y°, ^ M = 

In other words one has to check validity of the conditions (Y"_ n )_„ = (Y_„)_ u , (^l/n)^ = (\I/„)u, 
{y-u)-v = (y- v )-u, (i^Dv = i^PDu- Moreover, it is enough to take j = 0, because for the other 
j it is really true. The conditions required are equivalent the following five equations: 



(35) 
(36) 
(37) 
(38) 
(39) 



X„ — X- 



Jv-f 



Z V A/ = l + (z-l) 



-i ■ 



X v — X- 



Jv-f 



ze 



A/ 



EX 
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The second and fourth of them are the consequences of the boundary condition z\ = z~ x and 
the others have already been proved to be consistent under the boundary condition. 

In order to find the continuum limit as h — > rewrite the boundary condition z\ = z~ x in 
a more explicit form 

f- 2 . v = r 1 + f- j\ v . (40) 

Hence f u = f + hD x f + 0{h 2 ) and f v = f + hD y f + 0(h 2 ) then setting h — ► one gets 
f- 2 = f- 1 + f - f 1 . Then it follows from © that D x D y (f + f- 1 ) = 0. Consequently, 
/ = —f-i+a(x)+b(y). Remove functions a(x) and b(y) by the dilatation / = f+a(x)/2+b(y)/2. 
As a result one gets the boundary condition searched f~ l = —f. 

Notice that when the additional constraint f(u, v, k) = f(u + v, k) is imposed on the chain 
(JTJ), which reduces it to 1+1 dimensional discrete Toda chain, the boundary condition (|4*UJ) is 
reduced to the form f(u, —1) = —f(u, 0) found earlier by Yu.Suris in [TU] (see also [TT]). 

Consider now the other half of the chain which is located on the left half- line k < k$. 
Formulae (fTBj) allow one to exclude from the Lax pair the functions y° and and then rewrite 
the shifted variables yzl and in the form 

yZl = y-'e^-f- 1 - i^ 1 , C 1 = r 1 ^ 1 -^ - iy" 1 . (41) 

Introduce additional dynamical variables ^J^ 1 := ipzl and Y^ 1 := y~ x . Then the previous 
equation can be transformed as follows 

VZl = y^ert-r 1 - ^-i = ^-leAT 1 -/" 1 _ ly-i. ( 42 ) 

Below we will use the following analog of Lemma 1. 

Lemma 3. Solutions of the equations (JT§jl and (J2UJ) can be chosen to satisfy the constraint 

R UV S = -ly. (43) 

To prove the lemma one has to express H and X through R and S and substitute them 
into 

Lemma 4.. The shifts FjJ, YrJ, ip" 1 , ^Z x u of the variables F" 1 and ^>~ l are 

linearly expressed through the finite number of elements of the enlarged dynamical set 

{y-\Y-\4)- l ,y- l ,4)- 2 ,y- 2 ,...}: 

Y-u = y~\ 

% Ruv 

Y ^ = z v c^ u '~ v ~ ^ u Y~^ gf„ _„—f_„ i 

K 1 = r\ 

= \ gf-u,-v~f-v — 1 _| \ g/-u,-u - / 1 y~ 1 , 

^— It,— 1! "J— ti,— t! 
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Lemma 4 is proved similarly to the lemma 2. 

Introduce infinite dimensional vectors-columns P and Q as 



/ ••• \ 

y~ 3 
y- 2 
y- 1 

) 



and infinite dimensional matrices 

/ 

A = 



Q 



I ■■■ \ 

U _1 J 



1 











— 2 f-3 
— u ■' 


1 











_ u J 


1 











1 






B = 



f" 3 -/-3 



-1 

t-2_f- 





-1 




















— V J 









v ••• 















Z- 




J it, — u 




( . 












\ 










1 














c = 




• • e? v 




-3 


1 




















x 'f- 2 1 











v • 









1 





/ 





o/it / 



-1 







f -i_f-i 

OJ u J 



D — ■■■ e^'~ r '' _ j 


V ••• o o o 

and two more matrices each with only two nonzero entries 



/ ■ 



v • 

Proposition 4 









R 



-v r 1 -r 1 



Ru 



\ 



( ■■■ 



\ ■■■ 



The following system of equations 



Z—nC 







,-1 , 

J u, — v J 






-^e J 



P-u = AP, P- v = BP + bQ, 
Q v = CQ, Q U = DQ + dP 



is consistent if and only if the function / = f(u,v, k) solves the semi-infinite lattice 

-v-i 



AS 



AS' 



z 



for 



oo < j < —1. 



(52) 
(53) 



4 Finite chains 



Close the semi-infinite chain (JHBJ), (pMj) by imposing an additional (degenerate) boundary con- 
dition e? N+1 = 0. The obtained chain coincides with (Jlj). One can close also the Lax pair by 
setting y N+1 = 0, ^ N+1 = 0. The Lax pair found can be represented in the matrix form 



P_ u = AP + aQ, P_ t , = BP, 
Q V = CQ + cP, Q u = DQ. 

with the following matrix coefficients of dimension (N + 2) x (N + 2) : 



.4 



( z u 
1 

e f -"- f 1 



B 





V o 
/ o 








V o 

( z v 




















pN—1 rN — 2 

oJ — u J 









fN _j-JV-1 
_ u J 



\ 





1 J 



J-v-S 








fl _fl 














rN—1 _ rN—1 







-1 

fN _ fN 

3 J — v J 



\ 



C 




1 








D 



eti-f 1 


V o o o 

( 1 

e f »- f -1 
r 















JV fiV-1 









V o 

















f N-l_ fN—l 
dJu J 



\ 





1 J 

\ 





-1 

fN_ fN 
dJ u J 



(54) 
(55) 



(56) 



(57) 



(5f 



(59) 
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The matrices a and b are also of the dimension (N + 2) x (N + 2). Their entries are all zero 
except the first two entries of the first columns: aoo = X- U)V Z- U , a w = coo = H_ UtV z v , 

c io — H v . 

If one imposes the condition (JT5|l on both ends of the chain, one gets the chain 



Af- 1 



AP 



Af 



4 ' 



IT' 



for 0<j<iV-l, 



(see chain (JBJ) in Introduction). Its Lax pair is given by 



P_ u = AP + AaQ, 
Q V = CQ + X^cP, 

where A is the spectral parameter and 



P- v = BP + bQ, 
Q U = DQ + dP, 



I 



A 












V o 





1 

f 1 - f 











1 








B 



( o 








V o 







e f-v~f _l 

pJ — V J 



c 





























1 

ofv / 











1 












1 

riV-2_ fN-3 








-1 

rN — 2 rN — 2 

dJ —v J 


















1 

fN-l_ fN—2 

DJ —u J 










1 
1 









r 



1 

V-2_yJV-3 






rN—1 _ rN—1 











1 

N-l _jN-2 



\ 







o J 











AT f N-l_ f N-l 







z_ v e J 








1 
1 



\ 







J 



(60) 



(61) 
(62) 



(63) 



(64) 



(65) 
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D 



( 







V o 



Ju-f 










/ o 




V 





V o 











D J U J 















1 

R N 

N-l rN-1 









fN-2_ f N-2 
oju J 














fN-1 _ rN-1 

D J U J 















f N-X_ f N-l 
2 g J it , — v J —v 



\ 





H v 




o / v o o 
\ / o ... 



u, — v * — v 



d 



) 



o 
o 











\ 





\ 




1 

■ S N 



N-l rN-1 
Ju 



where R N and S are defined through equations 



ON 



1L Z g J U,—V •> —v 



R 



N 



s 



N eJ ~fu 



qN Q 



r)N N f N -f N 

R z_ u e J - u J - u '0 



A" 



yN-1 



tdN cN 

n, uv o 



z N -1 



(66) 



(67) 



(68) 



(69) 
(70) 

(71) 



Explain how the spectral parameter has been introduced. Evidently each of the functions 
H and X is defined by the equations ()23|) . U P to the constant multiplier. Due to the 
constraint (|22|) imposed above, the only constant remains free, and it is taken as the spectral 
parameter. 
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